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ABSTRACT 

The polarization position angles (PA) of pulsar radio emission occupy a dis- 
tribution that can be much wider than what is expected from the average linear 
polarization and the off-pulse instrumental noise. Contrary to our limited under- 
standing of the emission mechanism, the excess dispersion in PA implies that pul- 
sar PAs vary in a random fashion. An eigenvalue analysis of the measured Stokes 
parameters is developed to determine the origin of the excess PA dispersion. The 
analysis is applied to sensitive, well-calibrated, polarization observations of PSR 
B1929+10 and PSR B2020+28. The analysis clarifies the origin of polarization 
fluctuations in the emission and reveals that the excess PA dispersion is caused 
by the isotropic inflation of the data point cluster formed from the measured 
Stokes parameters. The inflation of the cluster is not consistent with random 
fluctuations in PA, as might be expected from random changes in the orientation 
of the magnetic field lines in the emission region or from stochastic Faraday rota- 
tion in either the pulsar magnetosphere or the interstellar medium. The inflation 
of the cluster, and thus the excess PA dispersion, is attributed to randomly po- 
larized radiation in the received pulsar signal. The analysis also indicates that 
orthogonal polarization modes (OPM) occur where the radio emission is heavily 
modulated. In fact, OPM may only occur where the modulation index exceeds 
some critical value, (3 C ~ 0.3. 

Subject headings: Methods: data analysis - polarization - pulsars: general - 
pulsars: individual (PSR B1929+10, PSR B2020+28) 

1. INTRODUCTION 

The linear polarization of pulsar radio emission is generally attributed to charged parti- 
cles streaming along dipolar magnetic field lines above the pulsar's polar cap (Radhakrishnan 
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k, Cooke 1969). Given the pulsar's enormous magnetic field strength 10 12 gauss) and 
the stability of its magnetic field structure, as implied by the long term stability of pulsar 
average profiles (Helfand, Manchester, & Taylor 1975), one might expect the polarization 
position angle (PA) at a particular pulse longitude to be fixed rigidly on the plane of the 
sky. However, PA histograms constructed from single-pulse polarization observations are 
much broader than what is expected from the average linear polarization and the off-pulse 
instrumental noise (Manchester, Taylor, & Huguenin 1975, hereafter MTH; Stinebring et al. 
1984, hereafter SCRWB; McKinnon & Stinebring 1998, hereafter MSI; Karastergiou et al. 
2002), suggesting that pulsar PAs also vary in a random fashion. 

This excess dispersion in PA may be intrinsic to the pulsar, possibly arising from an 
analog of stochastic Faraday rotation in the pulsar magnetosphere or from randomly polar- 
ized radiation (RPR) in the emission. Since the orientation of a pulsar's linear polarization 
vector is thought to be aligned with the magnetic field lines in the emission region, the excess 
PA dispersion could be due to fluctuations in field line orientation. However, this explana- 
tion seems highly unlikely because it requires the orientation of a 10 12 gauss magnetic field 
to fluctuate on a timescale comparable with a pulsar's rotation period (P ~ Is). Stochastic 
Faraday rotation in the interstellar medium (ISM) also seems an unlikely explanation for 
the origin of the excess PA dispersion because the rotation measure (RM) of the ISM does 
not fluctuate on the timescale of a rotation period. Furthermore, RM fluctuations should 
cause PA variations to be similar everywhere within a pulsar's pulse, but the PA disper- 
sion appears to vary with pulse longitude. Whatever the cause of the excess dispersion, 
determining it requires a detailed understanding of what causes polarization fluctuations. 
To discriminate between mechanisms that create polarization fluctuations, analytical tools 
must be developed for a thorough interpretation of single-pulse polarization observations. 

The statistical model of radio polarimetry developed by McKinnon (2003b) provides a 
foundation for the detailed interpretation of polarization observations. The model accounts 
for the orthogonally polarized modes (OPM) in the emission and has been used to argue that 
OPM occur simultaneously, not separately (MSI; McKinnon & Stinebring 2000, hereafter 
MS2; McKinnon 2002). The model also accounts for deviations from mode orthogonality and 
has been used to show that the modal connecting bridge in PSR B2016+28 at 1404 MHz is a 
transition between modes of nonorthogonal polarization (McKinnon 2003a), thereby resolv- 
ing a possible discrepancy in the interpretation of the pulsar's viewing geometry (SCRWB). 
Additionally, the model has been used to show that the large fluctuations in fractional circu- 
lar polarization observed in pulsar radio emission (SCRWB) may be attributed to the heavy 
modulation of the mode intensities and the small degree of circular polarization intrinsic to 
the modes (McKinnon 2002). 
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The objective of this paper is to derive and apply the analytical tools necessary to deter- 
mine the origin of the excess PA dispersion in pulsar radio emission. An eigenvalue analysis 
of radio polarization measurements that is based upon the three-dimensional statistics of 
radio polarimetry (McKinnon 2003b) is presented in § 2. In the analysis, eigenvalues of 
polarization measurements are derived for different polarization models, such as fixed po- 
larization, RPR, random fluctuations in PA, OPM, and nonorthogonal polarization modes. 
The eigenvalue analysis is applied to sensitive, well-calibrated polarization observations of 
the nearby pulsars PSR B1929+10 and PSR B2020+28 in § 3. The analysis indicates that 
the excess PA dispersion is caused by RPR in the received pulsar signal. In § 4, the nature 
of polarization fluctuations in pulsar radio emission is clarified, the origin and consequences 
of RPR are discussed, and additional applications of the eigenvalue analysis are proposed. 
The conclusions of the analysis are summarized in § 5. 



2. Eigenvalue Analysis of Radio Polarimetry 

The Stokes parameters Q, U, and V completely describe the polarization of electromag- 
netic radiation. Fortunately for radio astronomers, all three Stokes parameters of a radio 
signal can be measured simultaneously, so the polarization state of the radiation can be com- 
pletely determined with a single measurement (Radhakrishanan 1999; McKinnon 2003b). 
When multiple measurements of the Stokes parameters are plotted in a three-dimensional 
Cartesian coordinate system, or Poincare space in what follows, the shape of the resulting 
data point cluster generally resembles an ellipsoid. The location of the ellipsoid's centroid 
is given by the average values of the Stokes parameters. The dimensions and orientation 
of the ellipsoid are determined by polarization fluctuations intrinsic to the radio source and 
by the instrumental noise. It follows that one may determine what causes the polarization 
fluctuations provided the dimensions and orientation of the ellipsoid can be measured. 

The dimensions and orientation of the polarization ellipsoid can be measured with an 
eigenvalue analysis of the Stokes parameters. The covariances of the three Stokes parameters 
form a 3x3 data matrix. The eigenvalues of this covariance matrix are measures of the ellip- 
soid's three dimensions. The eigenvectors associated with the eigenvalues are the three axes, 
or the orientation, of the ellipsoid. The dimensions of the ellipsoid are largely unaffected by 
whether the source polarization occurs in Q, U, V, or any combination of the three Stokes 
parameters. Therefore, the eigenvalue analysis and the conclusions drawn from it are gener- 
ally independent of whether the radio source is linearly-, circularly-, or elliptically-polarized. 
In the following examples, eigenvalues are derived for different types of mechanisms that 
create polarization fluctuations. 
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a) No Polarization 



b) Random Polarization 




c) Fixed Polarization 



.Q. 



d) Superposed OPM 



e) Disjoint OPM 



e- 



Fig. 1. — Shapes of Q-U-V data point clusters for different types of polarization. The shapes 
are shown as cross-sections of the cluster in the Q-V plane of Poincare space. The cluster 
is a sphere for no polarization (a), random polarization (b), and fixed polarization (c). The 
cluster resembles a prolate ellipsoid for superposed OPM (d). For disjoint OPM (e), the 
polarization data reside in two, diametrically-opposed ellipsoids. The signal in (c) and the 
OPM in (d) and (e) are linearly polarized. 
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2.1. Fixed Polarization 

Let us consider a radio source with constant linear polarization, /i, and assume that its 
polarization occurs in the Stokes parameter Q. When a radio telescope is pointed towards 
the source, the measured Stokes parameters are 

Q = /i + X N>Q , (1) 

U = X N)U , (2) 

V = X N>V , (3) 

where X N is a Gaussian random variable with zero mean and standard deviation, a N , that 
represents the instrumental noise. When multiple measurements of these Stokes parameters 
are plotted in Poincare space, the resulting cluster of data points is spherical in shape. The 
sphere has a radius that is characterized by the magnitude of the instrumental noise, cat, 
and its centroid is offset from the coordinate system origin by a distance \x (Fig. lc). 

The dimensions of the sphere can be determined through an eigenvalue analysis of the 
Stokes parameters. The matrix formed by the covariances of the Stokes parameters is 

/Cov(Q,Q) Cov(Q,U) Cov(Q,V)\ /l 0\ 

C= Cov(U,Q) Cov(U,U) Cov(U,V) = a 2 N 1 . (4) 
\Cov(V,Q) Cov(V,U) Cov(V,V)/ \0 l) 

Regardless of radio source polarization, the elements of the covariance matrix are always 
real and the matrix is always symmetric about its diagonal because the Stokes parameters 
are always real numbers (i.e., the covariance matrix is always Hermetian). In this particular 
example, the off-diagonal elements of C are equal to zero because the instrumental noise 
is independent between Stokes parameters (MSI; MS2). Since C is symmetric about its 
diagonal and its off-diagonal elements are equal to zero, the matrix forms the eigenbasis 
of the data point cluster, and the elements on the matrix diagonal are its eigenvalues. The 
square root of each eigenvalue is equal to the magnitude of the instrumental noise, a at, which 
is proportional to the size of the sphere's radius. The actual radius of the sphere is about 
3a N because the fluctuations in instrumental noise are Gaussian and, thus, approximately 
99.73% of the Q-U-V data points reside within a sphere of this radius. Since the eigenvalues 
are equal to one another, the matrix is triply-degenerate, which means that the eigenvectors 
corresponding to the eigenvalues do not uniquely define the three axes of the data point 
cluster, as one would expect for a sphere. 
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The matrix formed by the second moments of the Stokes parameters is 

/ 1 + s 2 \ 
S = all 10 




(5) 



where s = fi/cr N is the signal-to-noise ratio in the polarization vector amplitude. The matrix 
S is also an eigenbasis. The principal eigenvector of S has an eigenvalue of <J%(1 + s 2 ). It is 
aligned with the polarization vector of the radio source, which in this example is the Stokes 
parameter Q. 

If the radio source emits RPR, the Q-U-V data point cluster retains its spherical shape, 
but is inflated by an amount equal to the random polarization fluctuations, a P (Fig. lb). 
The eigenvalues of the covariance matrix become Tn = r 2 2 = t 33 = a 2 N + a 2 p , and the sphere's 
radius is proportional to the square root of these eigenvalues. If the random polarization 
fluctuations are also Gaussian, the measured standard deviation in PA is 



Stochastic Faraday rotation and random fluctuations in the orientation of a pulsar's 
magnetic field are possible explanations for the excess PA dispersion in the emission. Both 
processes can be modeled as random fluctuations in PA. For a linear polarization vector with 
fixed amplitude and an orientation that randomly varies in the Q-U plane of Poincare space 
(i.e. in azimuth only), the measured Stokes parameters are 



where is a random variable that accounts for the vector's fluctuations in azimuth. If the 
fluctuations in vector azimuth are Gaussian with zero mean and a standard deviation of a^, 
the measured polarization occurs primarily in Stokes Q, and the eigenvalues of the covariance 
matrix are 




(6) 



2.2. Random Fluctuations in Position Angle 



Q = /icos0 + X NjQ , 
U = /isin0 + A NjU , 
V = Xn,Vi 



(7) 
(8) 
(9) 



ni = ^-[l-exp(-aJ)] 2 + 4 
7-22 = y[l-exp(-2aj)] + ^ 



(10) 



(11) 
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t 33 = cr 2 N . (12) 

The shape of the Q-U-V data point cluster formed by Gaussian fluctuations in PA is a 
banana- shaped, arcing ellipsoid. The centroid of the arcing ellipsoid resides in the Q-U 
plane at a distance \x from the origin of Poincare space. When the PA fluctuations are small 
{(Ttj, <C 1), Tn — T33, and the cluster is a prolate ellipsoid (i.e., an ellipsoid that is rotationally 
symmetric about its major axis). The ellipsoid's major axis is aligned with the Stokes 
parameter U because the polarization fluctuations are largest in U (i.e., r 22 > Tu > r 33 ). 
Therefore, the major axis of the ellipsoid is perpendicular to the orientation of the radiation's 
polarization vector. Since polarization position angle is related to polarization vector azimuth 
by ip = 0/2, the standard deviation in the measured PA is 

** = \{s- a + o}) 1/a - (13) 



2.3. Orthogonal Modes of Polarization 

In the statistical model of pulsar polarization proposed by McKinnon & Stinebring 
(MSI; MS2), the observed polarization is determined by the simultaneous interaction of 
two, orthogonally polarized modes with randomly varying polarization amplitudes. Since 
the modes are orthogonal, their polarization vectors lie on the same diagonal in Poincare 
space, regardless of the degree of their linear or circular polarization. The statistical model 
designates the amplitudes of the mode polarizations by the random variables X\ and X 2 , 
which have means fi\ and yU 2 and standard deviations <j\ and <r 2 , respectively. Assuming 
for this particular example that the modes are linearly polarized and that their polarization 
vectors lie along the diagonal aligned with the Stokes parameter Q, the measured Stokes 
parameters are 

Q = X 1 -X 2 + X NiQ , (14) 
U = X NiU , (15) 

V = X N)V . (16) 

If Xi and X 2 are independent random variables, the covariance matrix computed from these 
Stokes parameters is 

/1+p 2 0\ 
C = o 2 N l 10, (17) 

V 1/ 

where p = (af + a 2 ) 1 ^ 2 /cr N is the magnitude of the intrinsic polarization fluctuations relative 
to the instrumental noise. The shape of the Q-U-V data point cluster formed by superposed 
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OPM resembles a prolate ellipsoid 1 (see Fig. Id). The centroid of the ellipsoid is offset from 
the origin of Poincare space by a distance p = \p\ — p 2 \. The dimensions of the ellipsoid 
are proportional to the square root of the eigenvalues. The size of the ellipsoid's major axis 
scales as ai = \fr xl = cr N (l + p 2 ) 1 / 2 . The dimensions of the ellipsoid's minor axes are scaled 
by the instrumental noise, a 2 = 03 = <Jn- The orientations of the ellipsoid's minor axes are 
not uniqely determined because the two minor eigenvalues are equal to each other. Unlike 
the case of random PA fluctuations considered above, the major axis of the OPM ellipsoid 
is parallel to the radiation's polarization vector. Thus, the OPM model predicts that the 
principal eigenvectors determined from the covariances and second moments of the Stokes 
parameters will be aligned. 

The axial ratio of the ellipsoid can be related to the modulation index of the emission. 
For example, if the modes are completely polarized, the polarization fluctuations are given 
by p = j3s, where (3 and s are the modulation index and signal to noise ratio, respectively, 
of the total intensity (McKinnon 2002). The ellipsoid axial ratio is 

a 1 /a 3 = (l+p 2 ) 1 / 2 = (l+/?V) 1 / 2 , (18) 

which becomes larger as the emission becomes more heavily modulated and as the signal to 
noise ratio of the observation improves. At the other extreme when the signal to noise ratio 
is poor or the polarization does not fluctuate (i.e. p — f3 — 0), the axial ratio approaches 
unity, and the cluster is a spheroid. 

The detailed shape and dimensions of an OPM Q-U-V cluster can be affected by the 
character of the fluctuations in X\ and X 2 . Simple numerical simulations of equations 14 
through 16 indicate that the cluster resembles a prolate ellipsoid whenever o\ ~ a 2 , regardless 
of the character of the fluctuations. The simulations also show that Gaussian fluctuations 
in Xi and X 2 produce a Q-U-V cluster that is a prolate ellipsoid, even if <j\ differs from 
o 2 . When X\ and X 2 are exponential random variables with <j\ ^> a 2 , the ellipsoid is still 
rotationally symmetric about its major axis, but the ellipsoid is stretched along its major axis 
in the direction of the dominant mode, X ± . This is to be expected because the polarization 
fluctuations would be largely determined by the dominant mode. The dimensions of the 
cluster's minor axes extend about ±3a 2 and ±3a3 from the cluster centroid because the size 
of the minor axes are set by the Gaussian instrumental noise. The size of the cluster's major 
axis extends about ±3ai from the centroid if the mode polarizations are Gaussian random 
variables and can be larger if they are exponential random variables. 



1 Even if the mode polarization amplitudes are not independent, the cluster is still a prolate ellipsoid but 
with a principal eigenvalue of Tn = o\ + erf — 2ri 2 ericr2 + <J%, where r i2 is the correlation coefficient of the 
polarization amplitudes (MSI). 
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If OPM are disjoint (i.e. occur separately) instead of superposed, the Q-U-V data points 
will reside in two ellipsoids that are diametrically opposed in Poincare space (Fig. le). The 
size of the ellipsoids' minor axes are set by the instrumental noise. The size of the ellipsoids' 
major axes are determined by the instrumental noise and the polarization fluctuations of the 
modes. 



2.4. Nonorthogonal Polarization Modes 

Nonorthogonal polarization modes do not occur along the same diagonal in Poincare 
space. However, the two vectors representing the modes still define a plane, and the orienta- 
tion of this plane depends upon the elliptical polarization of the modes. Here, and without 
loss of generality, the modes are assumed to be linearly polarized. If the deviation from or- 
thogonality is given by the angle 9, incorporating the nonorthogonality of the modes in the 
model described above (McKinnon 2003a) gives the following expressions for the measured 
Stokes parameters. 

Q = Xi — X 2 cos 9 + X Nj q, (19) 

U = X 2 sin# + X NiU , (20) 

V = X N , V . (21) 

If the mode polarizations are independent (i.e., Cov(X 1 ,X 2 ) = 0), the covariance matrix for 
the nonorthogonal modes is 



C = 



' a\ + of cos 2 9 + a 2 N 
—a\ cos 9 sin 9 




>f cos 9 sin 9 



of sin 2 9 




(22) 



The eigenvalues of the covariance matrix in equation 22 are 



Til 



(o 2 + of + 2a 2 N ) + ^ + o| + 2o 2 of cos(2fl) 



(23) 



T22 



K 2 + 4 



2a 



NJ 



af + al + 2a\a\ cos(2#) 



T33 = 



(24) 
(25) 



The shape of the Q-U-V data point cluster formed by nonorthogonal polarization modes also 
resembles an arcing ellipsoid because the nonorthogonal mode breaks the symmetry of the 
cluster. 
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3. APPLICATION OF EIGENVALUE ANALYSIS 

The eigenvalue analysis developed in § 2 was applied to sensitive, well-calibrated, po- 
larization observations of PSR B2020+28 and PSR B1929+10 to determine the origin of 
the excess PA dispersion. The observations were made by SCRWB with the Arecibo radio 
telescope at a frequency of 1404 MHz. At the time the observations were made, the tele- 
scope's nominal system temperature was 40 K, and its forward gain was 8 K/Jy. Data from 
the same observations were used in the analyses documented in MSI, MS2, and McKinnon 
(2003a). SCRWB's superb polarization calibration is essential to the success of the eigenvalue 
analysis because improper or inadequate calibration will alter the dimensions of the Q-U-V 
data point cluster, and thus adversely influence the conclusions drawn from the analysis. To 
insure that the observed polarization fluctuations were intrinsic to the pulsar, and not due 
to slowly- varying signal amplification caused by diffractive interstellar scintillation (MSI), a 
subset of the overall time series where the running mean intensity did not vary was selected 
for further analysis. The time series thus selected contained 800 pulses for PSR B2020+28 
and 4000 pulses for PSR B1929+10. The pulse profiles and modulation indices computed 
from the selected time series are shown in Figure 2. 

The eigenvalue analysis was performed as follows. For both pulsars, a covariance matrix 
of the Stokes parameters Q, U, and V was computed at each pulse phase bin. The covariances 
were then used to calculate the eigenvalues and eigenvectors of the matrix. As illustrated by 
the examples in § 2, the dimensions of the Q-U-V ellipsoid at each phase bin were estimated 
by taking the square root of the eigenvalues. Examples of Q-U-V ellipsoids for each pulsar are 
shown in Figures 3 and 5. In each panel of both figures, every data point has been rotated 
by the PA at the relevant phase bin so that the major axis of the polarization ellipsoid 
occurs in the Q-V plane of Poincare space. The rotation eliminates projection effects in the 
data presentation. The bottom panels of Figure 6 show how the ellipsoid axial ratios vary 
across the pulse of each pulsar. The ratio of the ellipsoid's largest dimension to its smallest 
dimension, ai/a 3 , is shown by the solid line in the bottom panels of the figure. The ratio of 
the ellipsoid's two smaller dimensions, a 2 /a 3 , is shown by the dotted line in the same panels. 

Eigenvalues and eigenvectors were also computed from the second moments of the Stokes 
parameters to insure that this new data analysis method reproduced the results obtained 
from more traditional methods. As predicted by the eigenvalue analysis, the PAs of polar- 
ization vectors given by the second moment principal eigenvectors were very similar to the 
PAs computed in the traditional way, 

. 1 

w — - arctan 
2 

Furthermore, the vector colatitude given by the second moment principal eigenvector was 
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Fig. 2.— Pulse profile and intensity modulation indices of PSR B2020+28 and PSR 
B1929+10 at 1404 MHz. The average profiles of the two pulsars are shown in the top panels. 
Total intensity, linear polarization, and circular polarization are denoted by the solid, dashed, 
and dotted lines, respectively. The bottom panels show how the total intensity modulation 
index varies across the pulse of each pulsar. 
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very similar to the colatitude traditionally calculated from 

V 



e 



arccos 



(q 2 + u 2 + v 2 y/\ 



(27) 



3.1. PSR B2020+28 

Figure 3 shows four examples of Q-U-V data point clusters in PSR B2020+28 and 
provides the clue that explains the excess PA dispersion in the pulsar. Data recorded off the 
pulse are shown by the cluster at phase bin 10. The cluster has the spheroidal shape and 
size expected for instrumental noise. The cluster's values of a±, a2, and are equal to each 
other and are equal to the noise calculated off the pulse, gn — 0.1 Jy. These values suggest 
that the cluster's three dimensions are equal and, consequently, that the cluster is spheroidal 
in shape. All the Q-U-V data points at bin 10 reside within a sphere of radius, r = 3a^. 
But contrary to most examples given in § 2, the smallest dimension of every other cluster 
shown in Figure 3, as indicated by the values of a 3 , is much greater than what is expected 
from the instrumental noise. 

At phase bin 74, the cluster is a spheroid that is three times larger than what is expected 
from the off-pulse noise. The fact that the cluster at bin 74 is a spheroid, and not an arcing 
ellipsoid, suggests that the excess PA dispersion at this pulse location is not caused by 
random fluctuations in PA (e.g. stochastic Faraday rotation or fluctuations in the orientation 
of magnetic field lines). Since the data point cluster is not a prolate ellipsoid, the data do 
not support the assumption of MSI and MS2 that both polarization modes occur at bin 74 
of PSR B2020+28. The large size of the spheroid, as well as the large values of a 3 at bins 61 
and 91, suggests that the excess PA dispersion is caused by a mechanism that isotropically 
inflates the clusters. The only mechanism considered in § 2 that can isotropically inflate the 
cluster is RPR in the received signal. 

Another intriguing and noteworthy feature of the data at bin 74 is that the fluctuations 
in total intensity and polarization are Gaussian-like (see Fig. 4 and Figs. 4 & 7 of MSI). This 
is in stark contrast to intensity fluctuations in other pulsars, which can be lognormal, power 
law, and chi-squared or gamma (e.g., Cairns, Johnston, & Das 2003; Cairns et al. 2003; MSI; 
Cordes 1976a, 1976b). The emission's extremely low modulation index near the center of 
the pulse is also indicative of Gaussian intensity fluctuations. Since total intensity must be 
non-negative, Gaussian intensity fluctuations require the mean intensity, /if, to be at least a 
factor of three to five times greater than the intensity standard deviation, 07. Therefore, the 
modulation index of Gaussian fluctuations must satisfy (5 — uj j < 0.2 — 0.3. Interestingly, 
at phase bins 65 to 85 where the modulation index is (3 < 0.3, the Q-U-V data point clusters 
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Fig. 3. — Q-U-V data point clusters at different locations in the pulse of PSR B2020+28. 
The pulse phase bin is denoted in the upper left corner of each panel. The quantities a±, a,2, 
and CZ3 in the bottom left corner of each panel are the square roots of the cluster eigenvalues, 
in units of Jy, and are proportional to the size of the cluster's three dimensions. For display 
purposes, the data points have been rotated so that the major axis of the cluster occurs in 
the Q-V plane of Poincare space. OPM occur at bins 61 and 91, but only one polarization 
mode occurs at bin 74. The overall size of each cluster is much larger than the off-pulse noise 
at bin 10, suggesting that randomly polarized radiation accompanies the radio emission. 
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also have spheroidal shapes, similar to that of bin 74. Thus, the polarization properties at 
bins 65 through 85 are similar to those at bin 74. 

The simplest explanation for the Q-U-V data point cluster at bin 74 is a combination 
of a fixed polarization vector, RPR, and instrumental noise. If this hypothesis is correct, 
the observed Gaussian fluctuations in polarization would be due to a combination of RPR 
and instrumental noise because the polarization vector is assumed to have a fixed amplitude. 
Furthermore, the distributions of PA, linear polarization, and circular polarization should 
all be broadened by the fluctuations due to RPR. To test this hypothesis, histograms of 
linear polarization, circular polarization, and PA constructed from the bin 74 polarization 
data were compared to the theoretical distributions expected for a fixed polarization vector 
in Gaussian noise (McKinnon 2002), where the noise is due to a combination of RPR and 
instrumental noise. 



UW = - exp 

71 



2a 2 



1 + 



|^cos(2^)exp(^cos 2 (2^) 



1 + erf 



/ fj, L cos (2jj) 
U V2 



x 



fy(x) 



exp 



1 



(// 



2a 2 



( XfJ, L 



aV27v 



exp 



(x - fj, Y y 
2a 2 



(28) 
(29) 

(30) 



The parameters used to construct the theoretical distributions were obtained from the ex- 
perimental data. The mean circular polarization was fiy — —0.32 Jy, the mean linear 
polarization was hl = (fJ>n + ^u) 1 ^ 2 = 1-H Jy and the standard deviation of the Gaussian 
noise was taken from the bin 74 data point cluster as a = 0.30 Jy. As shown in Figure 4, the 
theoretical distributions compare extremely well with the measured histograms. The com- 
parison of the PA distributions is much more favorable than the result obtained by MSI 2 
where identical data were used but the standard deviation of the Gaussian noise was as- 
sumed to be a = 0.1 Jy, instead of 0.3 Jy, in the calculation of the theoretical distribution. 
In summary, the favorable comparison between the theoretical and measured distributions 
in Figure 4 confirms a simple interpretation of the emission's polarization properties at bin 
74. The Gaussian fluctuations in polarization are due to RPR and instrumental noise. These 
fluctuations broaden the PA distribution as well as the distributions of linear and circular po- 
larization. Therefore, the excess dispersion in PA is due solely to RPR. Another component 
of the emission has a polarization vector with fixed amplitude and orientation. Consequently, 
only one polarization mode with constant polarization occurs at bin 74. 



2 The histograms in Figs. 7b and 7d of MSI utilize smaller bin sizes than those of Fig. 4. Therefore, the 
histograms appear to be slightly different, although the same data were used in the two figures. 
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Given this interpretation, one may also say something about the fluctuations in total 
intensity at bin 74. Since the mode's polarization is constant, the mode's intensity must also 
be constant because it is highly unlikely that the mode's intensity can randomly fluctuate 
while its polarization remains constant, particularly when the mode is completely polarized 
as assumed in MSI and MS2. Therefore, the observed fluctuations in total intensity at bin 
74 cannot be attributed to the polarization mode. They must be due to RPR. 

Gaussian fluctuations and random polarization are properties of instrumental noise. 
However, the inflation of the Q-U-V data point clusters cannot be attributed to on-pulse 
instrumental noise. At bin 74, the average total intensity is 2.7 Jy, or an antenna temperature 
of 21.6 K. This change in antenna temperature would increase the noise by no more than a 
factor of 1.5 (see eqn. (B4) of MS2), far short of the factor of three required to explain the 
size of the cluster. Furthermore, different locations within the pulse of PSR B2020+28 have 
similar total intensities, but the smallest dimensions of the clusters at these locations are very 
different (Fig. 6), suggesting that the mechanism that inflates the cluster is not instrumental 
in origin. Also, the inflation of the clusters cannot be attributed to the correlation of on- 
pulse instrumental noise between Stokes parameters. Correlated instrumental noise will 
stretch a cluster in one dimension, instead of inflating it, because correlated data points will 
tend to form a straight line in Poincare space. For example, in the extreme case when the 
instrumental noise is completely correlated between the Stokes parameters, all elements of 
the covariance matrix will be a 2 N . The matrix eigenvalues will be tu = 3<t^,T22 = t 33 = 0, 
which indicate that the data point cluster is one- dimensional. 

Superposed OPM definitely occur at bins 61 and 91 of PSR B2020+28 because the 
data point clusters at these locations are elongated along the line that passes through the 
cluster centroid and the Q-V origin (the diagonal in Poincare space). The simple fact that the 
cluster at each bin resembles a single, prolate ellipsoid and not two distinct ellipsoids suggests 
that OPM in this pulsar are superposed, not disjoint. At bin 61, the modes are primarily 
linearly polarized because the ellipsoid's major axis does not tilt much from V = 0. If the 
polarization data at bin 61 were plotted in a PA histogram, the modes would be shown to 
occur with nearly equal frequency (see Fig. 5 of MSI) because the ellipsoid centroid roughly 
coincides with the Q-V origin and the number of data points on either side of Q = is 
about the same. At bin 91, the modes have a significant degree of circular polarization, as 
indicated by the ellipsoid's tilt. One mode would occur much more frequently than the other 
in a PA histogram because most data points have Q > 0. RPR also accompanies the modes 
at these pulse locations because the two smaller ellipsoid dimensions indicated by a 2 and a 3 
are approximately equal but are much larger than the instrumental noise. 

Figure 6 indicates where OPM occur in the pulse of PSR B2020+28 and shows how the 
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Fig. 4. — Distributions of polarization position angle (a), linear polarization (b), and circular 
polarization (c) at pulse phase bin 74 of PSR B2020+28. The histograms were constructed 
from measured polarization values. The solid, continuous lines in each panel are the theo- 
retical distributions expected for a fixed polarization vector and Gaussian noise. 
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magnitude of RPR may vary across it. The RPR magnitude, as estimated from the value of 
a 3 for each cluster and the off-pulse instrumental noise, 

<J P = {al-alY'\ (31) 

is shown in Figure 6a. The magnitude of RPR is small in comparison to the total intensity 
at every location within the pulse. In Figure 6c, the large ratios of the ellipsoids' major 
and minor axes, 01/03, indicate where OPM occur, primarily in the peaks and wings of the 
pulse. The intensity modulation index is large (Fig. 2c) at the same locations. The ratio of 
the ellipsoids' smaller dimensions is 02/03 — 1 everywhere in the pulse, indicating that the 
ellipsoids are generally rotationally symmetric about their major axes as predicted by the 
OPM statistical model. If the small deviations from axial symmetry are real, they may be 
due to a slight nonorthogonality in the pulsar's polarization modes (McKinnon 2003a). The 
axial ratios and the modulation indices are small near the pulse center where it is possible 
that only one polarization mode with fixed polarization occurs. 



3.2. PSR B1929+10 

Figure 5 shows examples of Q-U-V data point clusters in PSR B1929+10. Again, the 
cluster representing the data recorded off the pulse has the spheroidal shape and size ex- 
pected for the instrumental noise. The magnitude of the noise is smaller than that for PSR 
B2020+28 primarily because a larger bandwidth was used for the PSR B1929+10 observa- 
tion. The smallest dimension of each ellipsoid is about a factor of two larger than expected 
from the instrumental noise, suggesting that the signal received from the pulsar contains 
RPR. 

The clusters at phase bins 90, 119, and 125 are highly elongated, single ellipsoids, 
indicating the presence of superposed polarization modes. Since the data at each bin do 
not reside in two distinct clusters, the observations do not support the hypothesis of disjoint 
polarization modes. The superposed modes are primarily linearly polarized because the 
major axis of each ellipsoid is nearly aligned with V = 0. Although both polarization modes 
occur at bin 90, as indicated by its highly elongated Q-U-V cluster, only one mode will appear 
in a PA histogram computed from these data because the vast majority of data points have 
Q > 0. This simply means that one polarization mode is much stronger than the other (i.e., 
/-ti 3> 1^2)- At phase bins 119 and 125, both polarization modes will appear in PA histograms 
because data points in these Q-U-V clusters reside on either side of Q = 0. In both cases, one 
mode will occur more frequently than the other because more data points have Q > than 
Q < 0. Figure 5 of McKinnon (2003a) shows the PA histogram for bin 119. The polarization 
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Fig. 5. — Q-U-V data point clusters at different locations in the pulse of PSR B1929+10. 
The pulse phase bin is denoted in the upper left corner of each panel, and the square roots 
of the cluster eigenvalues are listed in the upper right corner. The data points have been 
rotated so that the major axis of the cluster occurs in the Q-V plane of Poincare space. OPM 
occur at each bin where the cluster is highly elongated. The modes at bins 119 and 125 are 
not truly orthogonal because not all data points reside near the same diagonal. Similar to 
PSR B2020+28, the overall size of each cluster is larger than the off-pulse noise (bin 10). 
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modes at phase bins 119 and 125 are not precisely orthogonal because not all the data points 
in the Q-U-V clusters reside near a diagonal in Poincare space. These clusters are consistent 
with what is expected for superposed, nonorthogonal polarization modes (see § 2.4 above 
and McKinnon 2003a). Historically, cases of nonorthogonal polarization modes have been 
documented with PA histograms (Backer & Rankin 1980; SCRWB; McKinnon 2003a), which 
give no indication that the nonorthogonality may also occur in circular polarization. When 
the Q-U-V data are reported as polarization ellipsoids, the nonorthogonality in circular 
polarization is readily apparent. 

Figure 6b shows how the magnitude of RPR, as estimated from equation 31, varies 
across the pulse of PSR B1929+10. As with PSR B2020+28, the RPR in PSR B1929+10 
cannot be attributed to on-pulse instrumental noise because the total intensity of the pulsar 
cannot increase the noise by the factor required by the smallest dimensions of the polariza- 
tion ellipsoids. The magnitude of RPR in PSR B1929+10 is much less than that in PSR 
B2020+28. 

The shapes of the polarization ellipsoids in the pulse of PSR B1929+10 are consistent 
with the statistical model of superposed polarization modes. As shown in Figure 6d, the axial 
ratios of the ellipsoids are generally consistent with the prolate ellipsoidal shape predicted 
by the model. Near the center of the pulse, the ratios of the ellipsoids' minor axes are 
slightly greater than one, indicating that the ellipsoids are not perfectly symmetric about 
their major axes. The slight asymmetry is most likely due to the nonorthogonality of the 
polarization modes. At phase bins 95 through 125, both the ellipsoid axial ratio and the 
intensity modulation index (Fig. 2d) are relatively constant at ai/a 3 ~ 6.6 and j3 ~ 0.55, 
respectively, indicating that both polarization modes occur where the emission is heavily 
modulated. Throughout the pulse of PSR B1929+10, the modulation index is not consistent 
with Gaussian fluctuations in total intensity. 

SCRWB specifically commented on the very narrow distributions of PA and fractional 
circular polarization in PSR B1929+10. The PA distribution was the narrowest of the pulsars 
they studied, and they attributed the high degree of PA stability to an extremely homoge- 
neous emission region. The supplemental information provided to SCRWB's interpretation 
by the eigenvalue analysis is any fluctuations in PA and fractional circular polarization in 
excess of the instrumental noise may be attributed to RPR in the received signal. 
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Fig. 6. — Randomly polarized radiation (RPR) and the relative dimensions of the polar- 
ization ellipsoids in PSR B2020+28 and PSR B1929+10. The top panels show the total 
intensity (solid line) and the magnitude of RPR (dotted line) for each pulsar. The bottom 
panels show the relative dimensions of the polarization ellipsoids for each pulsar. The solid 
line is the ratio of the largest to the smallest ellipsoid dimension, and the dotted line is 
the ratio of the two smaller ellipsoid dimensions. The ratios indicate that the ellipsoids are 
generally prolate, a result that is consistent with the statistical model of superposed OPM. 
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4. DISCUSSION 

Three-dimensional measurements can assume a variety of shapes and sizes, and axial 
and eigenvalue ratios are often used to characterize these types of measurements (Fisher, 
Lewis, & Embleton 1987). Obviously, any model that hopes to describe the data must be 
able to reproduce the fundamental shape of the three-dimensional data set. The shapes and 
axial ratios of the data point clusters in PSR B2020+28 and PSR B1929+10 are consistent 
with the predictions of the statistical model of radio polarimetry (McKinnon 2003b). The 
model predicts that the clusters are prolate ellipsoids when superposed OPM occurs and 
are spheroids when the polarization is fixed. Clearly, the data are not consistent with some 
other cluster shape (e.g., dual ellipsoids or an oblate ellipsoid) that might be predicted by 
some other model (e.g., disjoint OPM). 

In summarizing their comprehensive polarization observations, SCRWB favored a modal 
broadening mechanism as an explanation for the broad PA histograms they observed. How- 
ever, they also allowed that the histograms could be explained by the superposition of OPM 
and RPR provided that a technique could be developed to separate the random compo- 
nent from the modal emission. Their interpretation also required the random component 
to be small in comparison to the modal emission. The analysis presented in this paper is a 
technique that allows the component separation to be made. In the analysis, RPR appears 
as the isotropic inflation of the Q-U-V data point cluster, and OPM manifests itself as an 
elongation of the cluster along a diagonal in Poincare space. When the eigenvalue analysis 
is applied to the observational data, the random component is found to be smaller than the 
modal emission as indicated by the large axial ratios of the clusters (Fig. 6). Expressed 
mathematically, the superposition of OPM and RPR represents a general statistical model 
for a polarized signal received from a pulsar. 

Q = sin 9 cos o (Xx - X 2 ) + X PjQ + X N>Q , (32) 

U = sin O sin o (Xi -X 2 ) + X PiU + X NjU , (33) 

V = cos 9 (X 1 - X 2 ) + X PiV + X N)V . (34) 

Here, o and O are the colatitude and azimuth, respectively, of the OPM polarization vectors 
(McKinnon 2003b), and Xp is a random variable representing RPR. This is not to say that 
Xi, X 2 , and X P always occur at every pulse longitude of every pulsar. It is conceivable that 
any one, or even all, of these components could be absent from the received signal. 

Karastergiou, Johnston, & Kramer (2003) offer a slightly different interpretation of the 
polarization fluctuations in pulsar radio emission. They suggest that in addition to the OPM 
fluctuations that occur along a diagonal in Poincare space, the orientation of the diagonal 
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randomly varies. Their interpretation does not include RPR. The Karastergiou et al. model 
can be represented by equations (32)-(34) if the RPR terms, X P , are ignored and if the 
angles 9 Q and (p Q are allowed to be random variables. 

The phrases "randomly polarized radiation" , "unpolarized radiation" , and "randomiza- 
tion of position angle" describe very different polarization processes, but they are occasion- 
ally used interchangably in the pulsar literature, particularly in interpretative discussions of 
mode-separated pulse profiles, single-pulse polarization observations, and radiation depolar- 
ization mechanisms. The differences between these processes must be explicitly described 
to avoid confusion in the use of the phrases. The Stokes parameters Q, U, and V recorded 
in an observation of an unpolarized radio source are consistent with instrumental noise. In 
Poincare space, they occupy a sphere of radius r = 3ajy centered on the origin (Fig. la). 
While RPR is unpolarized on average, it is not the same as unpolarized radiation. The 
Stokes parameters recorded in an observation of a randomly polarized radio source will also 
appear as a sphere centered on the origin of Poincare space (Fig. lb), but the sphere's radius 
will be r = 3((Tp + o"^ r ) 1 / 2 , where a P accounts for the random polarization fluctuations intrin- 
sic to the source. The phrase "randomization of position angle" implies fluctuations in the 
PA of a linear polarization vector which are restricted to the Q-U plane of Poincare space. 
Polarization observations of a radio source with random PA fluctuations will appear as an 
arcing ellipsoid. The two largest dimensions of the ellipsoid will correspond to eigenvectors 
in the Q-U plane. The ellipsoid's centroid will be offset from the origin by a distance equal 
to the polarization of the source. 

The Gaussian intensity fluctuations and the random polarization observed at the center 
of the PSR B2020+28 pulse suggest a possible origin of RPR. Pulsar radio emission has 
been described as an incoherent addition of individually coherent shot pulses (Rickett 1975; 
Cordes 1976a, 1976b; Hankins et al. 2003). The addition can occur within the telescope 
receiving system, in the dispersion and scattering of the pulses as they propagate through the 
ISM, and at the pulsar if many independent emission events contribute to the instantaneous 
signal (Cordes 1976a, 1976b). From the Central Limit Theorem, the recorded signal will have 
Gaussian statistics regardless of the intensity distribution of the shot pulses, provided the 
number of pulses in the incoherent sum is very large. If the shot pulses also have no preferred 
polarization, the polarization of the recorded signal will be completely random. This scenario 
suggests that the intensity and polarization of the RPR component (Xj> in eqns. (32)- 
(34)) always follow Gaussian statistics. It is difficult to extract the statistics of X P in PSR 
B1929+10 and at other locations within PSR B2020+28 because of the large polarization 
fluctuations due to OPM. Cairns et al. (2003) have found evidence for a small Gaussian 
emission component in the Vela pulsar using theoretical models of its longitude-dependent 
intensity distributions. They also attribute the Gaussian component to the incoherent sum 
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of shot pulses. 

The eigenvalue analysis may be used to determine how pulsar radio emission depolarizes 
at high frequency (Manchester, Taylor, & Huguenin 1973; Morris, Graham, & Sieber 1981). 
The depolarization has been attributed to OPM, randomization of PA, a process intrinsic 
to the emission mechanism, or any combination of the three (MTH; SCRWB; Xilouris et al. 
1994; McKinnon 1997; MSI; Karastergiou et al. 2002). Since the observed Q-U-V clusters 
are not the arcing ellipsoids expected for fluctuations in PA, the analysis presented here 
suggests that randomization of PA does not occur. Therefore, randomization of PA is not a 
strong candidate for a depolarization mechanism. However, the analysis indicates that RPR 
occurs in the received pulsar signal, and RPR may be a viable depolarization mechanism if 
it is generated at emission. Applying the eigenvalue analysis to high frequency polarization 
observations of pulsars may determine if OPM, RPR, or both lead to the depolarization. 
OPM will be the depolarization mechanism if the high frequency Q-U-V clusters are highly 
elongated ellipsoids with dimensions of their minor axes set by the instrumental noise. Since 
the radiation is depolarized, the cluster centroid should be located near the origin of Poincare 
space. If RPR depolarizes the emission, the high frequency Q-U-V clusters will be spheroids 
with radii much larger than the instrumental noise and with centroids near the Poincare 
origin. Applying the eigenvalue analysis to multi-frequency, polarization observations of 
pulsars can also be used to measure the spectrum of RPR and the longitudinal extent of 
OPM, both of which will be useful in determining the depolarization mechanism. 

The presence of RPR in the received signal has other subtle, yet important, conse- 
quences for pulsar polarization observations and our interpretation of them. Perhaps the 
most important consequence of RPR is that it masquerades as on-pulse instrumental noise, 
thereby causing histograms of linear polarization, circular polarization, polarization position 
angle, fractional linear polarization, and fractional circular polarization to be much wider 
than what is expected from the off-pulse instrumental noise. McKinnon (2002) derived the 
statistics of fractional polarization in an attempt to explain the large dispersion in frac- 
tional circular polarization reported by SCRWB. The excess dispersion was attributed to 
a combination of the signal-to-noise ratio of the observations, the heavy modulation of the 
mode intensities, and the small degree of circular polarization intrinsic to OPM. However, 
the derivation did not account for RPR, which increases the effective noise and thus de- 
creases the effective signal-to-noise ratio. The factor of two to three increase in effective 
noise indicated by the eigenvalue analysis of PSR B2020+28 and PSR B1929+10 is enough 
to attribute the observed dispersion in fractional circular polarization to the reduced signal- 
to-noise ratio caused by RPR. Similarly, Kramer et al. (2002) found that the distribution of 
fractional circular polarization in the Vela pulsar was much wider than what was expected 
from the off-pulse instrumental noise. RPR could account for the excess noise in the emission 
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from Vela and, as mentioned above, may be the Gaussian emission component described by 
Cairns et al. (2003). McKinnon & Stinebring (2000) attributed excess polarization fluctu- 
ations in PSR B2020+28 to on-pulse instrumental noise, but their noise model could not 
account for all the observed fluctuations. The RPR revealed by the eigenvalue analysis can 
easily account for the excess fluctuations. RPR also causes the linear polarization to be 
slightly overestimated. When reporting the results of pulsar polarization observations, it is 
common practice to subtract a bias offset caused by the off-pulse instrumental noise from 
the computed values of linear polarization. While the bias value of 



is relevant for data recorded off the pulse, the actual value may be much larger on the pulse 
where RPR can occur. 



The OPM statistical model requires fluctuations in mode polarization amplitudes to ob- 
serve the switching between modes. Similarly, moding will not be observed if the polarization 
amplitudes do not fluctuate. If the polarization amplitude fluctuates, one might reasonably 
expect the total intensity to fluctuate, also (MSI; MS2). Consequently, the model predicts 
that OPM should occur where the radiation is heavily modulated (§ 2.3). This is clearly the 
case for PSR B2020+28 and PSR B1929+10, as shown in Figures 2 and 6. The presence 
of OPM elongates the polarization ellipsoid, giving it a large axial ratio. For both pulsars, 
the axial ratio is large where the modulation index is also high. The axial ratio in PSR 
B1929+10 is relatively constant across much of its pulse, and its modulation index remains 
constant over a similar range in pulse longitude. In PSR B2020+28, both the axial ratio of 
the polarization ellipsoid and the intensity modulation index are large at the same locations 
within the pulse, primarily in the pulse peaks and wings. As mentioned in § 3.1, the intensity 
fluctuations are Gaussian and the modulation index is /3 < 0.3 near the center of the PSR 
B2020+28 pulse. The polarization data in this region are consistent with the occurrence 
of only one polarization mode. This observation raises the intriguing possibility that both 
polarization modes do not occur where intensity fluctuations are Gaussian. Stated another 
way, perhaps both polarization modes occur only where the modulation index exceeds some 
critical value, f3 c ~ 0.3. 

The longitudinal extent of OPM determined from the eigenvalue analysis (Fig. 6) is 
different from what has been determined in previous analyses. MS2 calculated the mode- 
separated pulse profiles of PSR B0525+21 and PSR B2020+28. The model used in their 
analysis neglected RPR, and their results indicate that both modes occur at all locations 
within the pulse of each pulsar. However, the eigenvalue analysis presented in this paper 
reveals the presence of RPR and suggests that only one mode occurs near the center of the 



(L) = a Ny ^j2 



(35) 
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PSR B2020+28 pulse. Using PA histograms, SCRWB found that OPM extended over a 
limited range of pulse longitudes. The elongated Q-U-V clusters may be a better indicator 
of OPM because, like the bin 90 data in PSR B1929+10, the PA histograms may only reveal 
the presence of one mode although both modes actually occur. Therefore, the longitudinal 
extent of OPM may be slightly greater than what is documented in SCRWB. 

RPR complicates the mode separation technique presented in MS2. The technique as- 
sumes that two superposed modes, each completely polarized, are the only components of the 
emission. The simplicity of this assumption allows the mode intensities and polarizations to 
be calculated from the average values of the measured Stokes parameters. The occurrence of 
RPR introduces another variable in the model (e.g. Eqns. 32 - 34), and the mode separation 
is then indeterminate unless some simplifying assumption regarding RPR can be made. The 
dominant effect of RPR upon the results presented in Figures 4 and 5 of MS2 will be to 
reduce the intensities and polarizations of the modes at all pulse longitudes. This effect may 
be small because the magnitude of RPR estimated in this analysis appears to be small in 
relation to the total intensity of the emission (Fig. 6). 

The results of single-pulse polarization observations are notoriously difficult to summa- 
rize and report, owing to the extreme variability of pulsar polarization and the large volume 
of data. The analysis developed here is a new method for presenting and interpreting the 
results of these types of observations. It is significantly different from previous data presen- 
tation methods. The results of some of the first single-pulse polarization observations were 
reported by plotting the raw Q-U-V data for a series of single pulses (e.g., Lyne, Smith, & 
Graham 1971). These seminal observations showed that pulsar polarization can be quite 
high and highly variable. Although the orientation of the polarization vector is somewhat 
difficult to visualize with this data presentation method, it readily illustrates the polarization 
associated with the total intensity of a limited number of pulses. MTH reported the results 
of their single-pulse polarization observations as multiple polarization ellipses plotted across 
the total intensity profiles of about a dozen, consecutive, individual pulses. The radiation's 
polarization state is immediately obvious with the MTH data presentation method. These 
observations were some of the first to reveal conclusively that OPM occur in the emission. 
SCRWB and Backer & Rankin (1980) presented the results of their single-pulse polarization 
observations as histograms of fractional circular polarization, fractional linear polarization, 
and PA. The histograms, which are particularly well-suited for documenting the occurrence 
of OPM, are statistical summaries that conveniently illustrate how the polarization fluc- 
tuates across the pulse. The polarization data recorded for a large number of pulses can 
be reported with this method. However, the histograms do not indicate the instantaneous 
orientation of the radiation's polarization vector because the circular and linear polarization 
data are reported separately. The data presentation method developed in this paper focuses 
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on a detailed investigation of a very large number of polarization samples at a single pulse 
longitude at the expense of both the total intensity information and a more global, instanta- 
neous view of the polarization at other pulse longitudes. The method immediately shows the 
orientation of the radiation's polarization vector, reveals the presence of OPM through the 
elongation of the data point cluster, and indicates the presence of RPR through the inflation 
of the cluster. Data presented in this manner can be readily interpreted in the context of 
other mechanisms that cause polarization fluctuations. 

The objective of most radio polarization observations is to determine where the tip of a 
polarization vector resides in Poincare space. This is done simply by averaging the individual 
Stokes parameters Q, U, and V. The tip of the vector is the centroid of the Q-U-V data point 
cluster, and it can be located more precisely by reducing the cluster size (i.e. the instrumental 
noise) through increased bandwidth and longer integration times. But as shown with the 
eigenvalue analysis and the statistical model of radio polarimetry (McKinnon 2003b), more 
polarization information is available to the radio astronomer. The unique ability to measure 
all Stokes parameters simultaneously at radio wavelengths also allows the astronomer to 
determine the cluster's size, shape, and orientation. The physical processes that create 
polarization fluctuations can be investigated by measuring all geometrical properties of the 
cluster. 



5. CONCLUSIONS 

In summary, an eigenvalue analysis has been developed for the interpretation of pulsar 
polarization data. The analysis is based upon a three-dimensional statistical model of radio 
polarimetry and allows the linear and circular polarization of pulsar radio emission to be 
interpreted together, instead of separately as has been done typically in the past. The model 
is generally consistent with polarization observations of PSR B1929+10 and PSR B2020+28. 
The analysis clarifies the origin of polarization fluctuations in pulsar radio emission. The 
observed excess dispersion in PA is caused by the isotropic inflation of the Q-U-V data 
point cluster. The cluster inflation is not consistent with random fluctuations in PA and is 
attributed to RPR in the received pulsar signal. The presence of RPR may contribute to 
the depolarization of the emission at high radio frequency. A comparison of the model with 
the observational data suggests that OPM are superposed. OPM occur where the emission 
is heavily modulated, and may only occur where the modulation index exceeds a critical 
value, j3 c ~ 0.3. The analysis also reveals deviations from mode orthogonality in the circular 
polarization of PSR B1929+10. 
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